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O ' Abstract 

(N' 

Harrison, Perkins and Scott have proposed simple charged lepton and neutrino mass matrices 



> 

o 






that lead to the tribimaximal mixing C/tbm- We consider in this work an extension of the mass 
matrices so that the leptonic mixing matrix becomes C/pMNS = VlUt^-buW , where V[ is a unitary 
matrix needed to diagonalize the charged lepton mass matrix and W measures the deviation of 
the neutrino mixing matrix from the bimaximal form. Hence, corrections to C/tbm arise from 

^: 

D , both charged lepton and neutrino sectors. Following our previous work to assume a Qin-Ma- 

^^ 

like parametrization Vqm for the charged lepton mixing matrix V[ in which the CP-odd phase is 

(N 

^ ■ approximately maximal, we study the phenomenological implications in two different scenarios: 

'sj" . V[ = VXyi and V[ = Vqm- We find that the latter is more preferable, though both scenarios are 

'nI" \ consistent with the data within "ia ranges. The predicted reactor neutrino mixing angle ^13 in 

f^ ' both scenarios is consistent with the recent T2K and MINOS data. The leptonic CP violation 

characterized by the Jarlskog invariant Jcp is generally of order 10~^. 
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I. INTRODUCTION 

The large values of the solar (^12) and atmospheric (^23) mixing angles may be telling us 
about some new symmetries of leptons not presented in the quark sector and may provide 
a clue to the nature of the quark-lepton physics beyond the standard model. If there exists 
such a flavor symmetry in Nature, the tribimaximal (TBM) jl| pattern for the neutrino 
mixing will be a good zeroth order approximation to reality : 

sin^6'i2 = -, sin^6'23 = -, sin6'i3 = 0. (1) 

For example, in a well-motivated extension of the standard model through the inclusion of A4 
discrete symmetry, the TBM pattern comes out in a natural way in the work of 2|]. Although 
such a flavor symmetry is realized in Nature leading to exact TBM, in general there may 
be some deviations from TBM. Recent data of the T2K J3| and MINOS jj] Collaborations 
and the analysis based on global fits |5|, |6| of neutrino oscillations enter into a new phase of 
precise measurements of the neutrino mixing angles and mass-squared differences, indicating 
that the TBM mixing for three flavors of leptons should be modified. In the weak eigenstate 
basis, the Yukawa interactions in both neutrino and charged lepton sectors and the charged 
gauge interaction can be written as 

-C = ^WM, (ulY + hmdn + -^W- hY^L + H.c. . (2) 

When diagonalizing the neutrino and charged lepton mass matrices UlAi^U* = 
dia.g{mi,m2,Tns), UlmiUR = dia,g{me,m^,mr), one can rotate the neutrino and charged 
lepton fields from the weak eigenstates to the mass eigenstates u^ — ?> U^ul, ^l{r) -^ 
^HR)^L{R)- Then we obtain the leptonic 3x3 unitary mixing matrix f/pMNS = UlU^, from 
the charged current term in Eq. ([2]). In the standard parametrization of the leptonic mixing 
matrix f/pMNS? it is expressed in terms of three mixing angles and three CP-odd phases (one 
for the Dirac neutrino and two for the Majorana neutrino) [7| 
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Pu , (3) 



\ ■'23''12 ~ '-23'-12'513e " ■='23'-12 ~ L-23''12''13e L-23'-13 / 

where Sij = sin 9ij and Qj = cos 9ij, and P^, = diag(e*'^^, e**^^, 1) is a diagonal phase matrix 
which contains two CP-violating Majorana phases, one (or a combination) of which can be 



in principle explored through the neutrinoless double beta (Oz/2/3) decay 8|. For the global 
fits of the available data from neutrino oscillation experiments, we quote two recent analyses: 
one by Gonzalez- Garcia et al. |5| 
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and the other given by Fogli et al. with new reactor neutrino fluxes |6|: 
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The analysis by Fogli et al. includes the T2K (Sf and MINOS U results. The T2K Collab- 
oration |3| has announced that the value of ^13 is non-zero at 90% C.L. with the ranges 



0.03 (0.04) < sin2 2^i3 < 0.28 (0.34) 



(8) 



or 



4.99° (5.77°) < ^13 < 15.97° (17.83°) 



(9) 



for 5cp = 0, sin^ 2^23 = 1 and the normal (inverted) neutrino mass hierarchy. The MINOS 
Collaboration found 



with a best flt of 



sin^2^i3 < 0.12 (0.20) , 



sin2 2^i3 = 0.04i;H^I (0.079+°:°^^) 



(10) 



(11) 



for 5cp = 0, sin^ 26*23 = 1 and the normal (inverted) neutrino mass hierarchy. The experi- 
mental result of non-zero \Uez\ = sin ^13 implies that the TBM pattern should be modifled. 
However, properties related to the leptonic CP violation remain completely unknown yet. 



The trimaximal neutrino mixing was first proposed by Cabibbo [9|^ (see also 10|) 
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with oj = e^'^'^l'^ being a complex cube-root of unity. This mixing matrix has maximal CP 
violation with the Jarlskog invariant \Jcp\ = l/(6-\/3)- However, this trimaximal mixing 
pattern has been ruled out by current experimental data on neutrino oscillations. In their 
original work, Harrison, Perkins and Scott (HPS) [l| proposed to consider the simple mass 
matrices 
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that can lead to the tribimaximal mixing, where a, x, y and z are real parameters,^ Mf = 
mim\ and M^ = A^^A^J,. The mass matrices are diagonalized by the trimaximal matrix 
Vc for charged lepton fields and the bimaximal matrix Ubm defined below for neutrino 
fields, that is, V^MfVc = diag(mg, m^,m^) and UI^M^Ubm = diag(mf,m2,m|). The 
combination of trimaximal and bimaximal matrices leads to the so-called TBM mixing 
matrix: 
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^ The matrix originally given by Cabibbo was in the form 




If one considers A^ discrete symmetry, it will have two subgroups, namely, Z2 and Z3. The trimaximal 
matrix given in Eq. (fT2|) is obtained under Z3. 
^ Different from the choice of HPS, the matrix element y in Eq. (fT3|) can be in general introduced as 

who pointed 
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complex: e.g., (M^)i3 = y and (A/^)3i = y* . This case has been considered by Xing 
out that the off-diagonal terms in f/eM will acquire a phase from the complex y. It has the interesting 
implication that a nonzero sin0i3 will result from the phase of y. However, the corresponding Jarlskog 
invariant is exactly zero and the absence of intrinsic CP violation makes this possibility less interesting. 



It is clear by now that the tribimaximal mixing is not consistent with the recent experimental 
data on the reactor mixing angle 6*13 because of the vanishing matrix element t/gs in f^TBM- 
In this work we consider an extension of the tribimaximal mixing by considering small 
perturbations to the mass matrices Mj and M^ which we will call M^ ^ and M'^ ^, respectively 
(see Eq. (IT^ below) so that U^ = UbmW is no longer in the bimaximal form and Ul = VcV[ 
deviates from the trimaximal structure, where V[ is the unitary matrix needed to diagonalize 
the matrix V^M^ ^V^. As a consequence, f/pMNS = UlU^ = V^ f/xBMW^ = f^TBM+ small 
perturbations. Hence, the corrections to the TBM pattern arise from both charged lepton 
and neutrino sectors. Inspired by the T2K and MINOS measurements of a sizable reactor 
angle ^13, there exist in the literature intensive studies of possible deviations from the exact 
TBM pattern. However, most of these investigations were focused on the modification of 
TBM arising from either the neutrino sector 12| or the charged lepton part (l3|, [ij], but not 



both simultaneously. 

The paper is organized as follows. In Sec. II, we set up the model by making a general 
extension to the charged lepton and neutrino mass matrices. Then in Sec. HI we study the 
phenomenological implications by considering two different scenarios for the charged lepton 
mixing matrix. Our conclusions are summarized in Sec. IV. 

II. A SIMPLE AND REALISTIC EXTENSION 



In order to discuss the deviation from the TBM mixing, let us consider a simple and 
general extension of the original proposal by HPS given in Eq. flT3|) . by taking into account 
perturbative effects on the mass matrices Mj and M^. The generalized mass matrices M' 
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and Ml ^ can be introduced as 
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^ TBM could be obtained in models with different discrete symmetries, such as S'3, A4, 5*4, A5. dihedral 
groups, • • • , etc. By considering higher order and radiative effects, the matrices in Eq. (|T5|) can be 
realized. For example, we have shown in Ref. jf 5| that these matrices can be obtained by introducing 
dimcnsion-5 operators to the Lagrangians. 



' ^ ~ m'rm'l 



where M C ^ and M^ ^ are defined as the hermitian square of the mass matrices Mr 

and M^ ^ = Ai'^^Ai'^, respectively, with the subscript / denoting charged fermion fields 



(charged leptons or quarks). Due to the hermiticity of Mr and M^ , the parameters 
ci'y 91,2,3,^^0, x',y', p are real, while b and Xi,2,3 are complex. The parameters (71,2,3, Xi,2,3 and 
p represent small perturbations. Note that the (11), (13), (22) elements {i.e., mix', mly' 
and TTig) in M^ ^ are assumed to contain any perturbative effects on the elements x, y, and 
z in M^, respectively. For simplicity, it is assumed that y' is real just as the other elements 



in Ml and the vanishing off-diagonal elements in M^ remain zeros in Ml 
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The parameters a and h are encoded in [l| as 
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where the subscript /, indicates a generation of charged fermion field, and fhf. represents a 
bare mass of fi, for example, fhf^ = rhe ^ rhf,^ = ffi^ <^ ffif^ = rhr for charged lepton fields. 
We first discuss the hermitian square of the neutrino mass matrix. Ml ^, in Eq. (IT5l) . It 
can be diagonalized by 
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with 



tan26' = — - 



(18) 



and 
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W= 1 Pu, (19) 

I -(cos^-sin^)/V2 (cos^ + sin^)/v^ , 

where the diagonal phase matrix P^, contains two additional phases, which can be absorbed 
into the neutrino mass eigenstate fields. For a small perturbation \p\ (^ l^^'l), the mixing 
parameter 6 can be expressed in terms of 



9 = 7r/4 + e with |e| < 1 . 



(20) 



W is then reduced to 
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The neutrino mass eigenvalues are obtained as 
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and their differences are given by 
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from which we have a relation Ar/igi — i^'^ii — ''^2(1 ~ ^')- I^ is well known that the sign 
of Am^i is positive due to the requirement of the Mikheyev-Smirnov-Wolfenstein resonance 
for solar neutrinos. The sign of Am|]^ depends on that of p/ sin 2e: Amli > for the normal 
mass spectrum and A?Ti|^ < for the inverted one. The quantities m\,m\^m\^9 (or e) are 
determined by the four parameters ml,x',y', p, while the Majorana phases in Eq. flTTl) are 
hidden in the squared mass eigenvalues. 

We next turn to the hermitian square of the mass matrix for charged fermions in Eq. ( ITSl) . 
This modified charged fermion mass matrix is no longer diagonalized by Vc 
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m"., rn%, respectively, and rjij is composed of the combinations 



a + b + b* 

corresponding to ml 

of 5'i,2,3 and Xi,2,3- To diagonalize V^M'^'^Vc = Vl (\i&g{ni^j^,rn?j^,rn'^j,^ V£ , we need an 
additional matrix V[ which can be, in general, parametrized in terms of three mixing angles 
and six phases: 
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where Sj = sin^j, q = cos^j and a diagonal phase matrix Pf = diag(e*^\ e*^^, e*^^) which 
can be rotated away by the phase redefinition of left-charged fermion fields. The charged 
fermion mixing matrix now reads Ul = VcVl- 

Finally, we arrive at the general expression for the leptonic mixing matrix 

t/pMNs = UlU, = vI^UtbmW . (27) 

A simple and general extension of the mass matrices given in Eq. ( 1T5|) thus leads to two 
possible sources of corrections to the tribimaximal mixing: Vl measures the deviation of the 
charged lepton mixing matrix from the trimaximal form and W characterizes the departure 
of the neutrino mixing from the bimaximal one. The charged lepton mass matrix in Eq. ( !T5|l 
or flMj) has 12 free parameters. Three of them are replaced by the phases ^1,2,3 in Eq- f l26|) 
which can be eliminated by a redefinition of the physical charged lepton fields. The remaining 
9 parameters can be expressed in terms of rrie, m^, m^, 6*1, 6*2, O^,., tti, 0^2, 03. /,From Eqs. (|2] 
and f l26p the mixing angles and phases can be expressed as 

n 1^231 ^ |'7l3| ^ |'7l2| I ^ 

fc'2 ^ -nr , fc'3 ^ -nr 5 «i = arg(?723), 



771;; 771^ mz^ 

a2 ~ - arg(?723) + arg(?7i3), "3 ^^ - [arg(77i3) - arg(?723)] + arg(77i2) , (28) 

with the condition m^^ ^ '722,'7ii- In the charged fermion sector, there is a qualitative 
feature that distinguishes the neutrino sector from the charged fermion one. The mass 
spectrum of the charged leptons exhibits a similar hierarchical pattern to that of the down- 
type quarks, unlike that of the up-type quarks which show a much stronger hierarchical 
pattern. For example, in terms of the Cabbibo angle A = sin6'c ~ \Vus\i the fermion masses 
scale as {me,m^) ^ {X^,X'^)mr, {md,ms) ~ (A^, A^)mb and {mu,mc) ~ {X^,X^)mt. This 



may lead to two implications: (i) the Cabibbo-Kobayashi-Maskawa (CKM) matrix [16[ is 
mainly governed by the down-type quark mixing matrix, and (ii) the charged lepton mixing 
matrix is similar to that of the down-type quark one. Therefore, we shall assume that (i) 
VcKM = ^L and V^ = 1, where V^ (V^) is associated with the diagonalization of the down- 
type (up-type) quark mass matrix and 1 is a 3 x 3 unit matrix, and (ii) the charged lepton 
mixing matrix ¥[ has the same structure as the CKM matrix, that is, V^ = Vckm or V^j^j^. 
Recently, we have proposed a simple ansatz for the charged lepton mixing matrix ¥[, 
namely, it has the Qin-Ma-like parametrization in which the CP-odd phase is approximately 
maximal [l3|. Armed with this ansatz, we notice that the 6 parameters 6*1, 6*2, ^3, ai, 02, ^3 
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in Vl are reduced to four independent ones f,h,\, S. It has the advantage that the TBM 
predictions of sin^6'23 = 1/2 and especially sin^ 6'i2 = 1/3 will not be spoiled and that 
a sizable reactor mixing angle ^13 and a large Dirac CP-odd phase are obtained in the 
mixing f/pMNS = ^l f^TBM- The Qin-Ma (QM) parametrization of the quark CKM matrix 
is a Wolfenstein-like parametrization and can be expanded in terms of the small parameter 

the QM one has the 



A 



17l |. However, unlike the original Wolfenstein parametrization 



advantage that its CP-odd phase 6 is manifested in the parametrization and is near maximal, 
i.e., 5 ~ 90°. This is crucial for a viable neutrino phenomenology. It should be stressed that 
one can also use any parametrization for the CKM matrix as a starting point. As shown 
in [l9|, one can adjust the phase differences in the diagonal phase matrix Py in Eq. (1261) in 
such a way that the prediction of sin^ 6*12 will not be considerably affected. 



For V]J = Vqm, the QM parametrization 13|, [iTJ] can be obtained from Eq. fl26l) by the 
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On the other hand, for Vl = Vqm the QM parametrization is obtained by the replacements 
sie*°i = (/ + he-'^)X^ , S2 = /lA^ , S3 = A , ^2 = , 03 = ^ : 
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where the superscript / denotes d (down- type quarks) or d. (charged leptons). From the 
global fits to the quark mixing matrix given by 20| we obtain 



/ = 0.7491--, h 



S09+°°^'^ 



A = 0.22545 ±0.00065, 5 
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(31) 



Because of the freedom of the phase redefinition for the quark fields, we have shown in 21 1 
that the QM parametrization is indeed equivalent to the Wolfenstein one in the quark sector. 
Finally, the leptonic mixing parameters (^23, 6*12, 6'i3, ^c-p) except Majorana phases can 
be expressed in terms of five parameters 9 (or e), 6,f,h,X, the last four being the QM 
parameters in the lepton sector. If we further assume that all the QM parameters except 5 
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have the same values in both the CKM and PMNS matrices, then only two free parameters 
left in the lepton mixing matrix are e and 6. If 6 is fixed to be the same as the CKM one, 
then there will be only one free parameter e in our calculation. In the next section, we 
shall study the dependence of the mixing angles sin^6'23, sin^^i2, sin ^13 and the Jarlskog 
invariant Jcp on 6 and e. 

To make our point clearer, let us summarize the reduction of the number of independent 
parameters in this work. In the leptonic sector, we start with 16 free parameters (12 from 
the charged lepton mass matrix M^ ^ and 4 from the neutrino mass matrix M^ ^) as shown 
in Eq. f lT^ . Among the 12 parameters from M^ ^, three phases can be rotated away by 
the redefinition of the charged lepton fields. The remaining 9 parameters correspond to 
three charged lepton masses [me^^^r) and six angles in the charged lepton mixing matrix 
Vl as shown in Eq. f l26|) . while the 4 parameters from M^ ^ correspond to three neutrino 
masses (^^1^2,3) plus one angle {6 or e) in the neutrino mixing matrix Ui, as shown in Eq. ( !T7|) 
or ( !2T]) . With our ansatz for Vl discussed before, the 6 angles in ¥[ are reduced to four 
QM parameters (/, h, A, 6). Thus, the number of parameters finally becomes five (/, h, A, 5 
plus 6 (or e)), except for the six lepton masses. Under the further assumption of the QM 
parameters /, h, A having the same values in both the CKM and PMNS matrices, these five 
parameters are reduced to only two ones S and e. 

III. NEUTRINO PHENOMENOLOGY 

We now proceed to discuss the low energy neutrino phenomenology with the neutrino 
mixing matrix f/^, (see Eq. f lT7|l ) characterized by the mixing angle 6 or the small parameter 
e and the charged lepton mixing matrix Ul = VcVr in which Vl is assumed to have the 



similar expression as the QM parametrization 13|, ll7| given by Vqj^ or Vqm (see Eq. ( 129|) 
and Eq. fl30|l . respectively). The lepton mixing matrix thus has the form 

„ J VqmUtbmW for Vl^ = Vqm, 

[ V^^UtbmW for Vl = Vqm- 

Therefore, the corrections to the TBM matrix within our framework arise from the charged 
lepton mixing matrix Vl characterized by the parameters f,h,X,6 and the matrix W speci- 
fied by the parameter e whose size is strongly constrained by the recent T2K data. Indeed, 
the parameters A, /, h and 6 in the lepton sector are a priori not necessarily the same 
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as that in the quark sector. Hereafter, we shall use the central values in Eq. (13T!) of the 
parameters (A, /, h) for our numerical calculations. 
In the following we consider both cases: 



(i)V, 



''-Vc 



QM 



With the help of Eqs. flT^ and fl29|) . the leptonic mixing matrix corrected by the replace- 
ments Vc —^ Ul = VcV[ = VcVqm ^^*i U^{tx /A) — ?> Uy{ii /A + e), can be written, up to order 
of A'^ and e^, as 
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Note that f/pMNs ^^^^ contains five independent parameters (A, h, /, 5 and e).^ By rephasing 
the lepton and neutrino fields e — )• ee*"\ /i — )■ /le*^^, r — ?■ re*'^^ and z/2 — )■ z/2 e*^"^""^-*, the 
PMNS matrix is recast to 
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P. 
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where Uaj is an element of the PMNS matrix with a = e, fi,T corresponding to the lepton 
flavors and j = 1, 2, 3 to the light neutrino mass eigenstates. In Eq. fl34p the phases defined 
as tti = arg(f/ei), 0^2 = arg(f/e2), as = arg(f/e3), (3i = arg(t/^3) and (32 = arg(t/^3) have the 



Our previous work 13[ corresponds to case (i) with e = 
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expressions: 

ai = tan~^ 



^2 



tan ^ 



A{ V3(e2 _ 2) sin 5 + 6e cos 5 - QheX^} 
V3(2 - e2)(2 - A2 - /iA3) + v^(e2 - 2)A cos 6 - 6eA sin 5^ 

A sin (5 
1 + A cos 5 - ^ + /iA3 



3^ 



_i / A{2y3e sin (5 + 3(2 - e^) cos 5 - 3h{2 - e^)X^} 

Oi'i — tan ■= 

\3A(e2 - 2) sin 5 - 2v^e(2 - X'^ - \cos5 - h\ 

_i ^ 3(2 - 6^)(2 - A^ - 2/A^) - 6/i(2 - 6^)A^ cos^ - 4V36A(2 + hX) sm6 ^ 
~ ^^^ \ 2v^e(2 - A2 + 2/A2) + 4y3eA(2 + hX) cos 6 - 6h{2 - e^)X^ sin 6 

_^ f3{2-e^){l + fX^) + 3hX\2-e^)cos6 + 2V3eX\h-2fX)sm6\ 

Do — tan ;= — I . (35) 

y 2V2e(/A2-l) + 2y3eA2(/i + 2/A)cos(5-3/iA2(2-e2)sin(5 ' 

^Prom Eq. ( IMl) . the neutrino mixing parameters can be displayed as 

\TT |2 \TT |2 

,2/1 I'^e2| .Oq I^^mSI 



•2/1 ^ ez -2/1 



sin 6*13 = \Uez\ , 5cp = ai - 03 . (36) 

It follows from Eqs. ( 1M|l and (136|) that the solar neutrino mixing angle 612 can be approxi- 
mated, up to order A^ and e^, as 

2 ^ 1 2e^ 2A / ^ e sin (5 e^ cos 5 
sm 6'i2 — - H — - — h ^r I cos H 7= — h 



3 9 3V V3 3 

A^ /I 2e sin 2(5 -2 



+ ^(2h-'-^ + j(2h-7cosi)Y (37) 

This indicates that the deviation from sin2 6'i2 = 1/3 becomes small when cos 5 approaches 
to zero and the magnitude of e is less than A. Since it is the first column of V[ that 
makes the major contribution to sin^^i2, this explains why we need a phase of order 90° 
for the element {y[)2i'- When | sin(5| ~ 1, the present data of the solar mixing angle can be 
accommodated even for a large |e| (but less than A). The behavior of sin^ 612 as a function of 
5 is plotted in Fig. [1] where the horizontal dashed lines denote the upper and lower bounds 
of the experimental data in 3o" ranges. The allowed regions for 5 (in radian) lie in the ranges 
of 1.45 < (5 < 2.17 and 4.17 < 5 < 4.91 , recalling that the QM phase is 5qm = 1.56 . 
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Likewise, the atmospheric neutrino mixing angle ^23 comes out as 

• 2/1 1 ^^ f ■ ^ ecofi5\ 

sm^.3-2-7l(--^-^ 

- AM - + / + /icos5 + e ^ h ^-(1 - /- /i cos 5 - cos 25) 

v4 v3 3 



3 /sin 5, 



A'^e --={?> + Ah cos 6) -e 
V2v3 



Q Of 

■ COS 6 + h — 2h cos^ 5 



(38) 



6 

Fig. [1] shows a small deviation from the TBM atmospheric mixing angle with 6*23 < 45° for 

< |e| < A. Owing to the absence of corrections to the first order of A or e in Eq. (!38l) . the 

deviation from the maximal mixing of 623 comes mainly from the terms associated with A^ 

or eA. Especially, for sin5 ~ 1 we have the approximation sin^ ^23 — | ~ — -^ — A^(/ + |), 

which implies sin^ 6'23 < 1/2 for < |e| < A. We see from Fig. [1] that sin^ 6'23 lies in the 

ranges 0.43 < sin^ 6*23 < 0.45 for < |e| < 0.1 . 

The reactor mixing angle 613 now reads 

. 2^ 2eAsin5 2e^ 

sm 6'i3 = ■j= 1 — — (1 — Acoso) 

V 3 3 

2/1 2^ » "^ /sin (5 2he €cos6 



Evidently, sin 6*13 depends considerably on the parameters A and e. Thus, we have a non- 



vanishing 6'i3 with a central value of sin6'i3 = \/\/2 or ^13 = 9.2° for e = 13|. Note 
that the size of the unknown parameter e is constrained by the plot of sin 613 versus 6 in 
Fig. [1] where the horizontal dot-dashed lines represent the present T2K data for the normal 
neutrino mass hierarchy. For a negative value of e, the plot for sin ^13 versus 6 is flipped 
upside-down. Assuming p > 0, we see from Eq. ( 1231) that a positive (negative) value of e leads 
to a normal (inverted) neutrino mass spectrum. For example, we flnd -j= < sin6'i3 < 0.22 
(0.07 < sinei3 < ^) for 5 = 1.56 and e < 0.08 (e > -0.11) . 

Leptonic CP violation can be detected through the neutrino oscillations which are sen- 



sitive to the Dirac CP-phase Sep, but insensitive to the Majorana phases in f/pMNS 2M- ^^ 



follows from Eqs. 0351) and 0361) that the Dirac phase 6cp = ai — 03 has the expression 
_i , a/3A{-2 cos (5 + a cos 25 + A^ cos 5 + VSeX sin 25} 



JCP 



tan 



2 VSA sin 5{1 - A cos 5 - f } - 4e{l - A cos5 - fA^ cos2 5 - \%h - ^)} 



(40) 
where terms of order e^, A^, e^A^ have been neglected in both numerator and denominator. 
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Assuming p > 0, we show in Table [J the predictions for 5p and 9i^ as a function of e , where 
we have used the central values of Eq. fl3Tl) . 



TABLE I: 


Predictions of 5cp and 6*13 as a function of e in the case of VjJ - 


= Vqm- 


e 




6cp [deg.] 


^13 [deg.] 


-0.012 ~ 0.08 
-0.11 0.012 




-173.6 ~ -169 
184.6 ~ 186.4 


9.4 ~ 5.8 
14.4 ~ 9.4 



To see how the parameters are correlated with low energy CP violation measurable through 
neutrino oscillations, let us consider the leptonic CP violation parameter defined through the 



Jarlskog invariant Jcp = Im[f/ei^M2^e2^/!i] = | sin 2^12 sin 26'23 sin 2^i3 cos 6'i3 sin (5cp (23| 
which is expressed as 

A / . ^ 4 cos 5 



JcP = -^^--[s.n6 + e^-2e^s.n5] (41) 



9 



( sin 6{h + cos 5) — ev3(l + / — cos 25 + hcos5) — e^(/i + 4 cos 5) sin 5) . 



We see from the above equation that Jcp is strongly correlated with e and 6 for the fixed 
values of A, h and /. As long as e 7^ (associated with the neutrino part) or A 7^ 
(associated with the charged lepton part), Jcp has a non- vanishing value, indicating a signal 
of CP violation. Eq. ( 14T|) could be approximated as Jcp ~ "TT? ~ |sin(5. The behavior 
of Jcp is plotted in Fig. [1] as a function of 6. When sin 5 ^ 1, it is reduced to Jcp ~ 
-^ - f < -|(> -f ) for e > (e < 0). Assuming p > 0, we find -0.050 < Jcp < -0.037 
(-0.037 < Jcp < -0.017) for e < 0.08 (e > -0.11) and 5 = 1.56. 

(ii) Vl = Vqm 

The resulting leptonic mixing matrix in this case can be expressed, up to order of A^ and 
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FIG. 1: Mixing angles sin 623, sin 612, sin0i3 and the Jarlskog invariant Jcp as a function of the 
CP-odd phase 5 for case (i), where the sohd, dashed and dot-dashed curves are for e = 0.01, 0.05 
and 0.1, respectively, the horizontal dashed lines denote the upper and lower bounds in 3a set by 
Fogh et al. (see Eq. (0)) and use of A = 0.2254, / = 0.749, h = 0.309 has been made. In the plot 
of sin 013 versus (5, the horizontal dot-dashed lines represent the present data of T2K Collaboration 
for the normal neutrino mass hierarchy, see Eq. ([8]). 



e^, as 



U. 



(ii) 
PMNS 



UlU,{'KlA + e) = V^JJ^^u 



W 



I 



^TBM + 



./- 

2 \/ 3 



+ A 



v^ 2^6 Ve 2^2 

9 ■ 9 

\ V2 "^ 2V6 ^ V6 2V2 / 



V6 



^(1+/A^) \ 



V3 



v^ 



2e-*'^ + (l+/+/ie-")A e-'-*-{i+/+/ie-")A i{i-/-he-'*)A 



v 



V6 
(/+fee''^)A-2A^fe 
V6 



V3 

(/+/ie'^)A+feA^ 

V3 



i(/+fee''^)A 
V2 



/ 



+ Ae 



ie''^(l+/A^) , -e''^£+Ae+/e''^A^£ 
v^ 2^6 






(42) 



2v% 



i(/+fce-''^'-|)A _ 2e-''^'£+(|+/+fee-'^)A£ 2e-''»Xi+/+fee-''^)A _^ i(/+fee~''^-|)Ae 



V^ 



2V6 



V6 



2V2 



t(/+fee''^)A I (/+fee'^)Ae-2feA^£ 



V2 



2V6 



(/+fee'*)A+2A^fc I i{f+he'^)Xe 



V6 



2^2 



15 



0.65 t 

0.60 

0.55 L 

^' 0.50 

i» 0.45 

0.40 1 

0.35 

0.30 




0.25 



3 4 

6 [rad] 



isns 





FIG. 2: Same as Fig. [T] except for case (ii). 



Just as in case (i), the exact TBM is recovered when both e and A go to zero. With the help 
of Eqs. fl36|) and fH2l) . the solar neutrino mixing angle 6*12 can be approximated as 

^2 



sin^ 612 ^ ^ + ^ ^ ( cos 6 + 



+ 



1 
3 
A2 



2e2 2A 

1 2e sin 26 
2 



y (2/COS5 



V3 3 

esin(5 



e sin 5 e cos S 
{3 + 4cos^6) 



(1-2/) 



'COS 5 



(2/ + 7) 



(43) 



v^ ^" "^' ■ 3 

which leads to, as in case (i), a tiny deviation from sin^ 6*12 = 1/3 when cos 5 — ?■ and 
A > |e|. As expected, since the second column related to e in the matrix Eq. (142|) is zero, the 
solar mixing angle is not affected to the first order of e. Because of a minus sign in front of 
the A cos 5 term, which constitutes the major correction to sin 6*12, the plot of sin^6'i2 versus 
6 (see Fig. [2]) is turned upside-down, contrary to case (i). When sin5 ^ 1, the present data 
of the solar mixing angle are well accommodated even for a large |e| (but less than A). The 
allowed regions for 6 lie in the ranges of 1.0 < 5 < 1.7 and 4.5 < 6 < 5.3 . This indicates that 
when the CP-odd phase 6 is near maximal, the data of sin^ 6*12 can be easily accommodated 
in case (ii) but only marginally in case (i). Hence, the precise measurements of the solar 
mixing angle in future experiments will tell which scenario is more preferable. 
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/^From Eqs. (l36l) and (H2l) . the atmospheric neutrino mixing angle 623 comes out as 
.9/1 1 ^^ f ■ r ecos5,^ „ ^x\ 

- AM / - /icos5 - e ^ ^ (2sm^5 + / + hcos5) 

V4 \/3 3 ' 

+ ^1 ^sin5(3-2/-4/icos5) 



e 
3 L 



;i + 2V3-6/-4/icos5)cos5 + 8/isin2 5-4/il . (44) 



Fig. [2] shows a small deviation from the TBM atmospheric mixing angle with 623 > 45°, 

recalling that 623 < 45° in case (i). It is thus crucial to have precise measurements of the 

atmospheric mixing angle in the future to see whether 623 < 45° or 623 > 45° in order 

to test different scenarios. For sin 5 ~ 1 the deviation from the maximal mixing of 623 is 

approximated as sin^ 6^23 — | ~ ~^ + A^(/ — |), which leads to sin^ 6^23 > 1/2 for < |e| < A. 

The behavior of sin^ 6*23 is plotted in Fig. [2] as a function of 6. Likewise, the reactor mixing 

angle ^13 can be written as 

■ 2n 26^ 2eAsin5. ^ /tt . rx ^2 A 2\ 

sm2^i3 = — + (ecos5- V3sm5) + AM --e^j 

3 /sin 5(1 -2/) e cos 5(1 + 2/) 



We find 0.07 < sinei3 < ^ (^ < sinei3 < 0.22) for 5 = 1.56 and e < 0.12 (e > -0.07) . 
The Dirac phase 6cp has the expression 



6c p = tan 



-1 



A{v^(2 - e2)(l - /A2) sin 5 - 6e(l + /A^) cos<5} 
V3(2 - e2){2 - A2 + (1 - /A2)A cos 6} + 6eA(l + /A^) sin 6 



_ ^^^-1 ( A{2v/36(1 - f\') sin 5 + 3(2 - e'){l + f\') cos 5} \ 
\3A(e2 - 2)(1 + /A2) sin(5 + 2V3e(2 - A2 + (1 - /A2) cos5) J 

Assuming p > 0, we show in Table Ull the predictions for 6cp and 6*13 as a function of e, 
where we have focused on the central values of Eq. ( 1311) . 



The strength of CP violation Jcp can be expressed in a similar way to Eq. fj4T]) 

A^ / \ 

+ — (sin5(/i — cos 5) + ev3(l — f — cos 25 — hcos5) — e^{h — 4 cos 5) sin5 j , 
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TABI 


^E II: Predictions of 6cp and ^13 as a function of e in the case of Vl - 


= Vqm- 


e 


6cp [deg.] 


^13 [deg.] 


~ 0.08 
-0.11 ~0 


-7.1 ~ -5.2 
-15.6 ~ -7.1 


9.1 ~ 12.6 
4.1 ~ 9.1 



which can be approximated as Jcp ~ ~l~K + t ^in^. When sin 5 ?« 1, it is further reduced 
to Jcp ~ "TT^ + I < f (> |) for e > (e < 0). Assuming p > 0, we see from Fig. |2]that 
0.014 < Jcp < 0.037 (0.037 < Jcp < 0.05) for e < 0.12 (e > -0.07) and 5 = 1.56 . 



IV. CONCLUSION 

In their original work, Harrison, Perkins and Scott proposed simple charged lepton and 
neutrino mass matrices that lead to the tribimaximal mixing J/tem- In this paper we con- 
sidered a general extension of the mass matrices so that the lepton mixing matrix becomes 
t^PMNS = Vi^ UtbmW. Hence, corrections to the tribimaximal mixing arise from both charged 
lepton and neutrino sectors: the charged lepton mixing matrix ¥[ measures the deviation 
of from the trimaximal form and the W matrix characterizes the departure of the neutrino 
mixing from the bimaximal one. Following our previous work to assume a Qin-Ma-like 
parametrization Vqm for ¥[ in which the CP-odd phase is approximately maximal, we study 
the phenomenological implications in two different scenarios: ¥[ = Vq^ and ¥[ = Vqm- We 
found that both scenarios are consistent with the data within 3a ranges. Especially, the 
predicted central value of the reactor neutrino mixing angle ^13 = 9.2° is in good agreement 
with the recent T2K data. However, the data of sin^ 9i2 can be easily accommodated in the 
second scenario but only marginally in the first one. Hence, the precise measurements of the 
solar mixing angle in future experiments will test which scenario is more preferable. The 
leptonic CP violation characterized by the Jarlskog invariant Jcp is generally of order 10~^. 

Acknowledgments 

This work was supported in part by the National Science Council of R.O.C under 
Grants Numbers: NSC-97-2112-M-008-002-MY3, NSC-100-2112-M-001-009-MY3 and NSC- 



18 



99-281 1-M-OO 1-038. 



[1] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002) 
|arXiv:hep-p h/0202074 1 ; P. F. Harrison and W. G. Scott, Phys. Lett. B 535, 163 (2002) 
|arXiv:hep-ph/02032'09] . 



[2 
[3 

[4] 
[5] 

[6 

[7 
[8 
[9 

[lo: 



X. G. He, Y. Y. Keum and R. R. Volkas, JHEP 0604, 039 (2006) [a rXiv:hep-ph/0601001 |. 

K. Abe et al. [T2K Collaboration], Phys. Rev. Lett. 107, 041801 (2011) [arXiv:1106.2822] 

[hep-ex]]. 

P. Adamson et al. [MINOS Collaboration], |ar2£iYillQ8.Q015J [hep-ex]. 

M. C. Gonzalez-Garcia, M. Maltoni and J. Salvado, JHEP 1004, 056 (2010) 

|arXiv:1001.4"524V 3 [hep-ph]]. 

G. L. Fogh, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, Phys. Rev. D 84, 053007 

(2011) t arXiv:1106.6"028l [hep-ph]]. 

K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

J. Schechter and J. W. F. Valle, Phys. Rev. D 25, 2951 (1982). 

N. Cabibbo, Phys. Lett. B 72, 333 (1978). 

L. Wolfenstein, Phys. Rev. D 18, 958 (1978); P. F. Harrison and W. G. Scott, Phys. Lett. B 

333, 471 (1994); R. N. Mohapatra and S. Nussinov, Phys. Lett. B 346, 75 (1995). 



[11] Z. Z. Xing, Phys. Lett. B 533, 85 (2002) [arXiv:hep-ph/0204049| . 

[12] Y. Shimizu, M. Tanimoto and A. Watanabe, Prog. Theor. Phys. 126, 81 (2011) 
|arXiv:110 5.2929 [hep-ph]]; N. Qin and B. Q. Ma, Phys. Lett. B 702, 143 (2011) 
|arXiv:1106.3284 [hep-ph]]; Y. J. Zheng and B. Q. Ma. larXiv:1106.4040l [hep-ph]: E. Ma and 
D. Wegman, Phys. Rev. Lett. 107, 061803 (2011) [ arXiv: 1106.42691 [hep-ph]]: X. G. He and 
A. Zee, Phys. Rev. D 84, 053004 (2011) [arXiv: 1 106 .4359) [hep-ph] ] : T. Araki, Phys. Rev. D 84, 
037301 (2011) [arXiv: 1 106 .52111 [hep-ph] ] : S. Morisi, K. M. Patel and E. Peinado, Phys. Rev. D 
84, 053002 (2011) |arXiv:1107.0696l [hep-ph]]: W. Chao and Y. J. Zheng. [a?Xw: 1107.0738 1 [hep- 
ph]; S. Dev, S. Gupta and R. R. Gautam, Phys. Lett. B 704, 527 (2011) [arXiv:1107.1125l 
[hep-ph]]; R. d. A. Toorop, F. Feruglio and C. Hagedorn, Phys. Lett. B 703, 447 (2011) 
|arXiv:1107.3486l [hep-ph]]. 

[13] Y. H. Ahn, H. Y. Cheng and S. Oh. larXiv:1105.4460l [hep-ph] (unpublished). 

19 



[14] S. Dev, S. Gupta and R. R. Gautam, Phys. Lett. B 704, 527 (2011) |arXiv:1107.lT25] [hep- 
ph]]; P. S. Bhupal Dev, R. N. Mohapatra and M. Severson, Phys. Rev. D 84, 053005 (2011) 
|arXiv:1107.2"378] [hep-ph]]. 

[15] Y. H. Ahn, H. Y. Cheng and S. Oh, Phys. Rev. D 83, 076012 (2011) f arXiv:1102.0879] [hep-ph]]. 

[16] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. Theor. 
Phys. 49, 652 (1973). 

[17] N. Qin and B. Q. Ma, Phys. Rev. D 83, 033006 (2011) iarXiv:1101.4729l [hep-ph]]. 

[18] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[19] Y. Koide and H. Nishiura, Phys. Rev. D 79, 093005 (2009) [ arXiv:0811.2839] [hep-ph]]. 

[20] CKMfitter Group, J. Charles et al, Eur. Phys. J. C 41, 1 (2005) and updated results from 



http://ckmfitter.in2p3.fr 



[21] Y. H. Ahn, H. Y. Cheng, S. Oh, Phys. Lett. B 701, 614 (2011) [ arXiv: 1105. 0450) [hep-ph]]. 
[22] G. C. Branco, R. Gonzalez Felipe, F. R. Joaquim, I. Masina, M. N. Rebelo and C. A. Savoy, 

Phys. Rev. D 67, 073025 (2003) |arXiv:hep-ph/0211001| . 
[23] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); D. d. Wu, Phys. Rev. D 33, 860 (1986). 



20 



